Cooperative Model of Bacterial Sensing 
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Bacterial chemotaxis is controlled by the signalling of a cluster of receptors. A cooperative model 
is presented, in which coupling between neighbouring receptor dimers enhances the sensitivity with 
which stimuli can be detected, without diminishing the range of chemoeffector concentration over 
which chemotaxis can operate. Individual receptor dimers have two stable conformational states: 

on . 

one active, one inactive. Noise gives rise to a distribution between these states, with the probability 
influenced by ligand binding, and also by the conformational states of adjacent receptor dimers. The 
two-state model is solved, based on an equivalence with the Ising model in a randomly distributed 
magnetic field. The model has only two effective parameters, and unifies a number of experimental 
findings. According to the value of the parameter comparing coupling and noise, the signal can be 
arbitrarily sensitive to changes in the fraction of receptor dimers to which ligand is bound. The 
counteracting effect of a change of methylation level is mapped to an induced field in the Ising 
model. By returning the activity to the pre-stimulus level, this adapts the receptor cluster to a new 
ambient concentration of chemoeffector and ensures that a sensitive response can be maintained 
over a wide range of concentrations. 

Q j. PACS numbers: 87.10.+e, 87.22.-q, 05.20.-y 



I. INTRODUCTION 

_o 

One of the reasons why we think living things are special is that they show awareness of the environment: they 
respond sensitively to stimuli and can adapt to changes in the surroundings. Such biological complexity is displayed 
' even by bacteria, which, in order to survive, have to be aware of their precarious environment where various conditions, 
i Pf ■ such as nutrient and toxin levels, acidity and temperature, may change rapidly Q. In order to discover universal 
principles, applicable at many levels of biological complexity, by investigating a simple system, Adler revived studies 
t-H , on bacterial chemotaxis which had been intensively investigated a century ago . Recent genetic engineering methods 
^ ■ have made it a paradigmatic system of cellular signaling and adaptation. 

A bacterium such as Escherichia coli or Salmonella typhimurium swims smoothly by rotating a bundle of helical 
fiagella counterclockwise, but tumbles chaotically if the flagella rotate clockwise. When it moves towards a higher 
concentration of attractant, such as aspartate, it tumbles less frequently. When it moves towards a higher concentration 
of repellent, it tumbles more frequently. So the bacterium performs a biased random walk towards an attractant and 
away from a repellent. This phenomenon is called chemotaxis jj). It is mediated by receptors with extracytoplasmic 
ON . sensing domains, connected by transmembrane helices to signaling domains in the cytoplasm. The receptors, which 
are predominantly dimeric, cluster at one pole of the cell ||. There are several types of transmembrane receptors, 
which respond to different chemoeffectors but use the same signaling pathway |9,|J|, as shown in FIG. 1. Each 
r/2 ■ receptor dimer is joined to two CheA kinase molecules, via two CheW proteins, forming a 2:2:2 complex. CheA is 
autophosphorylated at a rate that is greatly enhanced by the receptor. The phosphate is then passed from CheA 
to a regulator protein, CheY. When phosphorylated CheY diffuses to a rotatory motor, the probability of clockwise 
rotation of the motor, and consequently the tumbling frequency of the bacterium, increases. The binding of a 
, chemoeffector ligand to a receptor dimer can trigger a transmembrane conformational change which regulates the 
autophosphorylation of CheA; attractant binding decreases the rate, while repellent binding increases it. In this way, 
an extra-cellular stimulus, i.e. ligand binding to a receptor, can modify the tumbling frequency. It is generally thought 
■ that there exist two stable conformational states of the receptor dimer: an "active" conformation which corresponds 
to a very high rate of CheA autophosphorylation; and an "inactive" conformation which corresponds to a lower rate. 

An important feature of chemotaxis is that the tumbling frequency does largely depends on recent changes of the 
concentration of chemoeffector. This is achieved through an adaptation process, which returns the activity of the 
system to the pre-stimulus level after a period of time. Adaptation is assured by a feedback loop which involves another 
regulator protein, CheB. Like CheY, CheB also receives a phosphate group from CheA. Phospho-CheB mediates a slow 
demethylation of the receptor, countering the action of CheR, which promotes methylation. Attractant binding also 
makes the receptor a better substrate for CheR. Since methylation enhances the autophosphorylation rate of CheA, 
the change in the rate of phospho-CheY production is gradually reversed by the feedback, and the the tumbling 
frequency ultimately returns to the pre-stimulus level. 
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Each sub-unit of a receptor dimer consists of two helices. It has been proposed that the transmembrane signaling 
involves a scissor- like or pivoting motion of the pair of subunits fflO^ , or a piston- like motion involving a conformation 
change within just one subunit {) 
the aspartate receptor, which reveal that attractant binding induces a displacement of one helix, down towards the 
cytoplasm, while the other three helices are not detectably perturbed p2|Jl^ ]. 

The chemotactic response is extraordinarily sensitive; as little as a single attractant molecule can trigger a detectable 
motor response p^ , ^q ,p|. Recently, Bray, Levin and Morton-Firth suggested that this sensitivity might be related 
to the clustering of receptors on the surface of the bacterium. Without discussing any underlying mechanism, they 
considered the possibility that the binding of a single molecule ligand affects the activity of a number of receptors, so 
that the response is augmented |lq| . In this article, we present a physical model of collective signalling in a cluster of 
receptors. We propose that the cluster responds as an entity, as a consequence of nearest-neighbor coupling between 
individual receptor dimers. In our model, the influence of one dimer on another depends only on its activity, and 
not on whether it is liganded. Noise causes each of the receptor dimers to fluctuate between active and inactive 
states. When a given receptor dimer binds a molecule ligand, the probability of it being active is altered. Owing 
to the coupling, the probabilities of activity of adjacent receptor dimers is also modified and this effect propagates 
throughout the cluster. Thereby, the response to a stimulus is amplified. Moreover, ligand binding is a rapid process, 
while which of the receptor dimers are liganded is random, thus the above effect is averaged out, and the overall signal 
is a statistical average quantity as a function of the fraction of liganded receptor dimers. The model can be cast as an 
elementary neural network and reduced to the Ising model. Thus the paradigmatic system of cellular signaling and 
adaptation is related to its counterpart in statistical mechanics. The model provides a simple, unifying framework to 
understand a large amount of experimental data. Sensitivity to very small changes of concentration, together with 
the ability to respond to gradients over a broad range of concentrations, can naturally be attained. The model might 
be applicable to a variety of cellular signalling processes which demand a combination of low threshold of response 
and wide dynamic range. 

The organisation of the rest of the article is as follows. In Section II, we analyse various experimental results and 
argue for the necessity of taking into account inter-dimer coupling, as well as noise. The neural network-like model 
is constructed in Section III; it is solved in Section IV, by reducing it to the Ising model with a randomly distributed 
magnetic field. Adaptation, and subsequent signaling of the adapted system, is specifically discussed in Section V. 
Section VI contains discussions and a summary. 
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II. COUPLING AND NOISE 



A. COUPLING 



The existence of coupling among receptor dimers is indicated by a number of experimental results. First, cooperation 
among receptors in signaling and adaptation is hinted at by the fact that most of the chemotactic receptors cluster into 
a patch, located at one pole of the cell [f|jj|. From the viewpoint of evolution, we might formulate a useful biological 
principle: An attribute that exists most probably confers advantages over possible alternatives, especially if the latter 
have some apparent merit. In the present case, if there were no cooperation among receptors, a uniform distribution 
over the surface would be optimal in efficiency for capturing molecules |17| |. Since, in practice, they are found to 
cluster together, there is most likely an advantage due to this feature. Therefore, coupling among the receptors might 
well play a role in signaling and adaptation. Coupling among receptor dimers can certainly improve the sensitivity. 
It can amplify the signal generated by a stimulus, as has been anticipated by some authors. 

Second, it has been found that signaling can occur through receptor dimers that have been genetically engineered so 
that one subunit lacks a signaling domain [ p|]l6 — |2Cf | - As mentioned above, a conformational change of only one subunit 
has been observed in the crystal structure P4 13| . If the two subunits have essential differences and only one of them 
is involved in the transmembrane signaling, then inter-dimer coupling is inevitable to explain the experimental results 
on truncated subunits (with 50% probability the truncation would have been made on the signaling subunit, and no 
signaling would occur if the dimers act independently). However, there is also the possibility that the binding of ligand 
to one subunit automatically suppresses binding to the other subunit; then the transmembrane signal can always be 
generated with either subunit. In this case, inter-dimer coupling is not essential to explain the above experiments. 
Coupling is, however, necessary in order to reconcile the fact that dimers with a truncated subunit are functional with 
the favored mechanism of methylation, which requires interactions between two subunits of the cytoplasmic domain 

Third, it has been proposed, based on experiments, that at least in certain cases, receptor methylation is related 
to dimer-dimer interactions, i.e. methyltransferase bound to one dimer can methylate other dimers |21| raj. Another 
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support for coupling is the remarkable mobility of the P2 domain of CheA, which provides the docking site for CheY 
and CheB. This mobility can serve to amplify the phosphorylation signal Finally, a recent in vitro experiment 
by Liu et al. showed that kinase activation by a soluble signaling domain construct involves the formation of a large 
complex, with about 14 receptor signaling domains per CheA p3|] . This appears to be a strong support for the 
coupling among receptor dimers in vivo. 

B. NOISE 

A proper consideration of noise is important for several reasons. Thermal noise is certainly a significant issue since, 
for biological molecules, the energy barriers between distinct conformational states are generally comparable to kT. 
Thus there is a high probability of stochastic transitions from one conformation to another. Indeed, experiments have 
indicated that there is considerable thermal motion in receptors. Disulphidc trapping studies of the galactose/glucose 
binding protein have revealed spontaneous, large amplitude thermal fluctuations of the protein backbone structure 
[5]. 

Moreover, noise can also provide benefits. In the absence of noise, nearest-neighbour coupling amongst receptor 
dimers would cause the activity to spread across the whole array and would inevitably make the response to different 
stimuli indistinguishable . If noise is taken into account, individual receptors flicker between active and inactive states 
and the overall signal is a statistical average, which naturally varies for different numbers of liganded receptor dimers. 

With the consideration of noise, it is worth emphasising that the signaling process should be understood within the 
context of dynamic equilibrium: When the concentration of chemoeffector is stable, the conformational state and the 
level of methylation of each receptor dimer fluctuates microscopically, but the mean activity of the system remains 
unchanged. This equilibrium is shifted when the concentration of chemoeffector is changed. 

III. THE MODEL 

We study the total signal of the cluster of receptors as a quasi-equilibrium property corresponding to a certain 
concentration of chemoeffector and a certain level of methylation. This approach is justified by the wide separation of 
time scales in this system: Ligand binding and protein conformational transitions occur within milliseconds. Changes 
in protein phosphorylation occur on a time scale of ~ 0.1s. The much slower adaptation process, associated with 
modulations of methylation level, is on a scale of minutes ||J2^] . 

The quasi-equilibrium state of a dynamics is determined by the minimum of a noisy "energy function" (a Lyapunov 
function). This "energy" is not necessarily the actual physical energy, since it may describe an effective dynamics that 
"renormalizes" the underlying chain of physical processes. Similarly, the noise may not be due only to the temperature, 
but in the present case it mostly is. Such a description, which reduces degrees of freedom, is especially effective when 
the underlying physical processes are complicated, or unclear in detail. A typical example is Hopfield's neural network 
model |^J,^5|. Here we adopt this approach for chemotactic signaling, but with a different interpretation and with 
the emphasis placed on the determination of the equilibrium activity as a function of the external stimuli. 

Consider a lattice of receptors, whose basic unit is the receptor dimer, or equivalently, the whole receptor-CheW- 
CheA dimer complex. Initially, we consider a system that has adapted to zero concentration of the chemoeffector, 
and investigate the response when the concentration changes so that a fraction c of receptor dimers are bound to 
chemoeffector molecule ligands. The state of a receptor is a function of the effect of ligand binding and the states of 
the neighboring receptor dimers. Characterizing the state of receptor dimer i by a variable V, (according to recent 
observation, it is the vertical position of one of four helices, but we are open to possible re-interpretation according 
to experimental findings), and the effect of the ligand binding to receptor dimer i by Hi, most generally we have 

V t = V i ({Vj fti },{H j }), (1) 

where {Hj} denotes the set of Hj for j = 1, 2, • • •, {V^} denotes the set of all Vj for j ^ i. The natural assumption 
is that Vi is affected only by Hi and the states of the nearest neighbours. Furthermore, for the two-state model, in 
which Vi has only two possible values, V° or V 1 , the McCulloch-Pitts threshold model J2^,^] is a natural assumption. 
Thereby 

Vi = V-OT TijVj +Si- Ut), with i>(x) = | £o I < o > (2) 
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where J7j is a threshold value, and describes the coupling among receptor dimers, which is assumed to be nonzero 
only for nearest neighbours. We adopt the convention that V° is the active conformation and V 1 is the inactive one. 
Then Hi > for attractant binding, which tends to inactivate receptors and depress the autophosphorylation rate 
of CheA, thereby decreasing the frequency of tumbling. Conversely, H i < for repellent binding. 

It is well known that if = Tji and Tu = 0, clearly valid in the present situation, the dynamics is determined 
by a Lyapunov function (or Hamiltonian) |2q| , 

H = - T H v i v i - 12 HiVi + 12 UiVi > ( 3 ) 

<ij> i i 

where < ij > represents pairs of nearest neighbours. Taking into account the noise, which induces a state distribution 
which is nearly a Boltzmann distribution |25| ] , the problem reduces to the statistical mechanics of a system with the 
above Hamiltonian. 

In the simplest interpretation, the noise is purely thermal, f3 = 1/kT, and Eq. (^|) may be identified as the effective 
physical energy. According to recent observation, V° and V 1 are the two stable positions of one of the four helices. 
Therefore Hi and T i3 Vj are forces due to ligand binding and coupling, respectively. The "Zeeman energy" dependent 
on ligand binding is due to the free energy exchange with bound ligand. Similarly, the coupling energy is due to free 
energy exchange with the cytoplasm or membrane, which mediates the effective coupling. 

Eq. (H) may be transformed to the "spin" representation by writing Si = 2(Vi — V°)/AV — 1, where AV^ = 
(V 1 - V®). Then 

U = - 12 J v S * S i - 12 BA + + (4) 

<ij> i 

where = TijAV /4, Bi = HiAV/2, and Eq is a constant given a distribution of {Bi}. Hi is a "Zeeman energy" 
due to an effective "magnetic field" independent of {Bi}, which determines the equilibrium configurations in the 
absence of {Bi}, i.e. without ligand binding. Without loss of generality, we set TCi = 0. Thus in the absence of 
{Bi}, and if the noise is sufficiently high, Si is equally distributed between 1 and —1, and the "magnetization" is zero. 
In other words, it is assumed that there is no energy difference between the active (Si — —1) and inactive (Si = 1) 
conformations for an isolated, unliganded receptor dimer i. The physics does not change if this difference is set to be 
nonzero. Ligand binding shifts the energy difference to 2Bi. 

We have now reduced the model to an Ising model. The activity of the array of receptors corresponds to the average 
magnetization of a lattice of spins, and ligand binding of a receptor dimer corresponds to a local magnetic field at a 
lattice site: Bi = B if receptor dimer i binds a chemoeffector ligand, otherwise Bi = 0. If the fraction of liganded 
receptor dimers is c, then the value of Bi is randomly distributed between B and with probability 

p(Bi) = c5(Bi - B) + (1 - c)8(Bi). (5) 

This Ising model in a field bimodally distributed between and B is simpler than the so-called "random-field Ising 
model" [p7|-j30||, in which the possible values of the field are symmetric with respect to zero, and nontrivial results arise 
due to the fluctuation of the fields. In our case, the average of the field is nonzero, so there is a long-range order simply 
as the result of the explicit symmetry breaking. For Eq. (|J), Bi = cB, where the overbar denotes the average over 
disordered configurations. The fluctuation of the random distribution is ABi — y c(l — c)B. Consider the formation 
of a domain of size L in the ferromagnetic ground state. According to the central limit theorem, the average Zeeman 
energy is ~ L d cB, much larger than its fluctuation, which is ~ L d l 2 \J c(\ — c)B. Therefore the energy gain is always 
positive, and the fluctuation of the field cannot destroy long-range order. 



IV. SOLUTIONS OF THE MODEL 



The two-state model, which has been reduced to the Ising model in a randomly distributed field, as described by 
Eqs. (Q) and can be solved by the mean-field method. One may obtain the result simply by considering that the 
average magnetization, m = < ij >, where < • • • > denotes the thermodynamic average, is determined by the local 
field B\ = Y]j Jijm + Bi with the random distribution. Alternatively one may first obtain the free energy using the 
replica method, then calculate the average magnetization pN,E9|. It is found that m is the root of the equation 



m = tanh((ivJm + /3Bi) (6) 
2c 2(1 - c) 



I + exp[-2(/3i/ Jm + (3B)} 1 + exp(-2/3^Jm) 



(7) 
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while the the noisy Lyapunov function is 



F= ^vJm 2 -^{c\n[2cosh(f3vJm + f3B))} + {l-c)\n[2cosh(f3vJm)}} + E . (8) 
2 p 

Here, has been assumed to have a single value J for nearest neighbors, v is the number of nearest neighbors, and 
(3 is a characterization of the noise. 

The relation between the chcmocffcctor concentration and the activity of the system is now reduced to the m versus 
c relation, determined by Eq. (Q), since the activity of the system, here defined as the fraction of receptor dimers in 
the active state is A = (1 — m)/2, and the pure response to the stimulus, i.e. the change of the activity, is AA — m/2. 
Although Eq. (fjj) may possibly have more than one solution, the one corresponding to the lowest F is what we need. 
Approximate analytical solutions may be found in limiting cases, 

P cB if 3 _> o 

l-2(l-c)exp(-2/3z/J) -2cexp[-2(/3z/J + /?.B)] if j3 -> oo and B > (9) 
-1 + 2(1 - c) cxp (-20vJ) + 2cexp -2[/3i/J + (3\B\) if j3 -► oo and 5 < 

In general, the solution can only be obtained numerically. It can be seen that there are actually only two effective 
parameters in this model; one is a = /3vJ, the other is 7 = f3B. Owing to symmetry, it suffices to give results for 
B > 0. Solutions for typical values of parameters are shown in FIG. 2. First we choose a — 0.1, 0.8, 1.2; then for 
each a, the dependence of m on c is determined for 7 = 0.01, 0.01, 0.1, 1, 10, 100. Note that for c = 0, i.e., the Ising 
model without a magnetic field, a = 1 is the critical value dividing the "paramagnetic" phase, where m(c = 0) = 0, 
and the "ferromagnetic" phase, where there is a "spontaneous magnetization", m(c = 0) =/= 0. 
The quantitative measure of sensitivity, denoted by S, is half of the slope at c = 0: 



S 



|c =° 

1 1 

l+cxp [-2(!3v.Jm a +fiB)\ l+oxp(-2 / 3i/Jm ) 
-1 4/3f J cxp( — 2/3^ Jm a ) 
1 _ [l+cxp(-20i/.7m o )J :i 



(10) 



where mo = m(c = 0). It is clear that S can be made arbitrarily large by choosing appropriate value of BvJ so that 
the denominator in the above expression is arbitrarily close to 0. For mo = 0, 
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C _ l+cxp(-2/3B) 2 ( y 

which is directly tuned by the difference between (3vJ and 1, which is the critical value of phase transition for c = 0. 
The case with mo ^ is less favored since the range of possible m for different c could be diminished, making it 
more difficult to distinguish between different stimuli. Moreover, the sign of mo would be determined by that of the 
previous B, conflicting the fact that the pre-stimulus level is fixed. 

For a given a, S also increases with 7 = /3B, but with an upper bound. The fact that dm/ 'dj — > when 7 — > 00 
indicates that, if ligand binding has a strong enough effect, the response is independent of the exact value of 7. This 
provides a sort of stability for the effect of ligand binding. 

Thus good sensitivity requires fine tuning of the coupling: the greater the sensitivity demanded by the bacterium, 
the more accurately a = (5vJ has to be controlled. But 7 = (3B may vary widely without considerably affecting the 
response. This is reasonable, since the temperature range suitable for bacterial survival is rather restricted and, for 
a given bacterium, v J is a structural property, which could be optimized during evolution. On the other hand, the 
effect of ligand binding, B, depends on the external stimulus, which may vary considerably. 

As an exercise, our model may be applied to the puzzling situation in which both attractants and repellents are 
present |]. In this case, 

p{B{) = CrSiB, - B r ) + CaSiB, - B a ) + (1 - Cr - C a )S(Bi), (12) 

where c r and c a are the concentrations of the repellents and attractants, respectively, and B r and B a are respectively 
the measures of the attractant and repellent binding. Obviously the activity is dependent on both c r and c a . 
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V. ADAPTATION 



Now we incorporate into this model the delayed adaptation due to the change of methylation level. This may 
be achieved through an induced "field" with an opposite sign to that associated with ligand binding, so that the 
"magnetization" returns toward the pre-stimulus level. This assumption for the additivity of the effect of ligand 
binding and that of the change of methylation level is supported by the finding that in a receptor there is a region 
which gathers, integrates, and interprets the multiple inputs transferred by the transmembrane signaling domain and 
the methylated side chains, then transmits an output signal to the kinase regulation machinery M . 

Two points should be made clear. First, the time scale of the change of the level of methylation of the whole system 
is much longer than the microscopic time scale, so the "magnetization" can still be obtained as the equilibrium 
property of the noisy Lyapunov function, which quasi-statically changes with the level of methylation. Second, since 
the ligand binding occurs on a time scale much shorter than the time needed for adaptation to be completely achieved, 
we cannot simply change the value of B, but must introduce another "field". We denote this "induced field" by {Mi}, 
with the distribution 



p(Mi) = C^Mj - M) + (1 - c m )5(M i ), (13) 

where c m is the fraction of the receptor dimers whose original level of methylation is modified. The sign of M is 
opposite to that of B. Thereby the net field is Di — Mi + Bi with the distribution 

p(Di) = cc m 5{Di -B- M) + c(l - c m )5(A - B) 

+ (1 - c)c m <5(A - M) + (1 - c)(l - c m )S(D t ). (14) 

The equilibrium state can be obtained by replacing Bi in Eq. (|j) with Di. Adaptation is taking place if c m and/or 
M vary slowly with time. This gives rise to a time-dependent "magnetization" , which may return to zeo. To get 
a quantitative impression, by adopting the high noise limit (3 — > 0, it may be estimated that when c m B + cM = 
0, adaptation is completed; the "magnetization" returns to zero. Here, we shall simply assume that a molecular 
mechanism exists which ensures that the state of zero "magnetization" is an attractor of the dynamics, so that 
adaptation is exact. A more precise study of the adaptation process will be reported in the future. 

Once the system has adapted, suppose that the concentration subsequently changes from c to c + Sc. One can 
obtain the new activity by substituting in Eq. ( |l4| ) c + Sc for c, and the values of c m and M at which the adaptation 
was completed. In general, what is most important is the change of fraction of liganded receptor dimers since the last 
adaptation. Under high temperature approximation, the result is Eq. (9) with c replaced by Sc. Moreover, it can be 
seen that if 5c is negative, i.e if chcmocffector is removed, the effect is similar to the addition of a chemoeffector whose 
"field" has opposite sign. Therefore the removal of attractant is equivalent to repellent binding, and vice versa. This 
has indeed been observed in experiments p2| . 



VI. SUMMARY AND DISCUSSIONS 



In this article, we analyse relevant experimental results and draw the conclusion that both inter-dimer coupling 
and noise are crucial in the mechanism of chemotactic signaling and adaptation. The ratio between their measures, 
a = (3vJ is one of the two effective parameters in the cooperative model we construct. A second parameter is the 
ratio between the measure of the effect of ligand binding and that of the noise, 7 = (3B. The essential features 
due to the balance of coupling and noise are well captured by the paradigmatic model of statistical mechanics, the 
Ising model. We made an attempt to map the underlying mechanism of collective effects in chemotactic signaling to 
the Ising model in a randomly distributed field, with the distribution reflecting the ligand occupancy. To complete 
the mapping, we adopted the basis of Hopfield's neural network model. The great difference between time scales of 
the various chemical and mechanical processes makes it feasible to obtain the signaling level as a quasi-equilibrium 
property of a noisy Lyapunov function. This Lyapunov function describes the dynamics "renormalizing" underlying 
complexity. 

Our model provides the following picture. An individual receptor dimer has two stable conformational states, an 
active one that corresponds to a high rate of CheA autophosphorylation, and an inactive one that corresponds to a 
low rate. Noise gives rise to a distribution between these states and the partition is influenced both by ligand binding 
and by the conformational states of the neighboring receptor dimers. In the simplest interpretation, the noise is purely 
thermal, (3 = 1/kT, 2B corresponds to the shift of the energy difference between active and inactive states induced 
by ligand binding, and J measures the effective coupling energy between neighboring receptor dimers. The activity 
of the receptor cluster is a statistical average quantity. A change in the fraction of liganded receptor dimers causes 
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the total activity to change from the pre-stimulus level. But the level of methylation also changes, on a slower time 
scale. This causes an effect opposite to that induced by ligand binding. Consequently, the total activity ultimately 
returns to the pre-stimulus level. 

The coupling between receptor dimers naturally provides the sensitivity to small stimuli observed in experiments. 
Additionally, the noise makes the response to different values of concentration changes distinctive. Sensitivity to small 
changes in the environment requires fine tuning of the parameter a, but 7 may vary without considerably affecting 
the response. The equivalence between the removal of attractant and the addition of repellent, or vice versa, has a 
natural explanation. 

Among problems for further investigation are the effects of finiteness of the number of receptor dimers, potential 
randomness in the coupling, and features that might be lost in the mean field solution. Correlation between Bi, or 
Mi, at different sites i is also a possibility and might have useful consequences. The finite-size effect and the "random 
field" due to the change of methylation level may destroy the "spontaneous magnetisation" that exists for a > 1, thus 
relaxing the constraint on the precision to which a must be specified to give high sensitivity. The mean-field solution 
is least accurate when c — ► \, since the fluctuation of the field is y/c{\ — c)B, which increases to the greatest asc-> \. 
Thus maybe the sensitivity is lower at moderate values of the occupancy c, than at the extremes c — > and c — > 1. 
However, this is not necessarily a limitation. The fractional occupancy c is related to the ambient concentration of 
ligand [L] by 

(15) 



[L] + K d ' 

where K$ is the dissociation constant. Thus 

\L\K d S[L] S[L] 

dc= ( [L ] + W W = C(1 " C) W (6) 

Given that the bacterium probably needs to detect a relative change in concentration, Ttjj we see that the greatest 

sensitivity to a change in occupancy is demanded when c — » or c — » 1, and the least when c — > \. 

It is well known that the the two-state threshold neural network model is equivalent to a model with continuous 
variables in the high gain limit p4| ], with the Lyapunov function 

H = - TiiViVi + E ^ r 9i\V)dV - E HtVi, (17) 
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with Ui = g t 1 (Vi) determined by 



<ij> 



where it, is interpreted as the soma potential, while C is the input capacitance of the cell membrane. When Eqs. fll'i 
and ( [l8|) are adopted for the network of chemoreceptor dimers, Vi is a variable characterizing the stable conformation, 
i.e. the (vertical) position of the mobile helix of the receptor dimer, and u, is the instantaneous position. Thus Eq. 
( |l8| ) could be the equation of motion describing the transient process of the movement of the mobile helix, in response 
to a force Hi generated by ligand binding, as well as forces due to couplings with the neighboring receptor dimers. 
Of course, whether modification of ( p~8| ) is necessary depends on future experimental results. 

According to this interpretation, by measuring the force generated by ligand binding, H, and the displacement of 
the mobile (signalling) helix AV = {V 1 - V°), one may obtain the parameter B = H ■ AV/2. Note that H ■ AV is 
the work done by the force H , consistent with the identification of 2B with the shift in energy difference caused by 
ligand binding. Similarly, AJ/AV — TijAV is the force generated by the conformational change of one of the nearest 
neighbours. To make a rough estimation, we take typical values a = 0.5, 7 = 5 and 1/(3 ~ ApN ■ nm (assuming that 
the noise is purely thermal). Then B « 20pN ■ nm, v,J « 2pN ■ nm. The measured displacement is 0.16nm Jl2| |. It 
is found that the force resulting from ligand binding is about 2b0pN and the force due to coupling between a pair of 
nearest neighbours is about WpN. These orders of magnitude are quite reasonable. 

Since the continuum model can be realized in electric circuits, more insights might be provided from the viewpoint 
of system control, where negative feedback has been well studied. On the other hand, the analogy with the neural 
network model is possibly more than a mathematical one. From the viewpoint of evolution, there are common features 
between bacterial sensing and sensing of higher animals. Perhaps a primitive or ancestral neural network works in 
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chemotaxis. Adler writes: "The basic elements that make behaviour possible in higher organisms are also present 
in a single bacterial cell; they are sensory receptors, a system that transmits and processes sensory information and 
effectors to produce movement. Whether the mechanisms of any of these elements in bacterium are similar to those 
in more complex organisms remains to be established" 0. Margulis thinks: " Thought and behaviour in people 
are rendered far less mysterious when we realize that choice and sensitivity are already exquisitely developed in the 
microbial cells that became our ancestors" . We hope our approach is a small step in addressing these issues. 
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FIG. 1. A schematic illustration of the chemotactic signaling pathway. 



FIG. 2. The solution of the two-state model: "magnetization" m as a function of the fraction c of liganded receptor dimers. 
Here we assume there was no ligand bound previously. The three figures are for three typical values of the parameter ct = /3v J: 
(a) a — 0.1; (b) a — 0.8; (c) a — 1.2. In each figure, different plots are for different values of the parameter 7 = f3B: 
0:7 = 0.01; + : 7 = 0.1; □ : 7 = 1; x : 7 = 10; A : 7 = 100. Note that the critical pont, which separates "ferromagnetic" and 
"paramagnetic" phases, is a = 1. 
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